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ABSTRACT 



We present the z = 0.3—0.7 cluster X-ray luminosity function (XLF) determined from 
the Southern SHARC (Serendipitous High-redshift Archival ROSAT Cluster) survey. 
^ \ Over the luminosity range L ~ (0.3 — 3) x 10 44 erg s _1 (0.5 - 2.0 keV) the XLF is 

in close agreement with that of the low redshift X-ray cluster population. This result 
greatly strengthens our previous claim of no evolution of the cluster population, at these 
luminosities, at a median redshift of z — 0.44. 
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1. Introduction 

Xray clusters of galaxies are efficient tracers of the 
mass in the Universe and can be studied out to large 
redshifts. Prior to the launch of the ROSAT X-ray 
telescope, the only high redshift (z > 0.3) X-ray se- 
lected cluster sample was that of the EINSTEIN Ex- 
tended Medium Sensitivity Survey (EMSS; Henry et 
al. 1992; Gioia & Luppino 1994). Henry et al. (1992) 
used the EMSS to show X-ray clusters evolved 'neg- 
atively' — the space density of high luminosity clus- 
ters being lower in the redshift range z = 0.30 — 0.60 
compared to z = 0.14 — 0.20. This result was in con- 
flict with popular models of cluster formation (e.g. 
Kaiser 1986) and prompted further significant theo- 
retical work (e.g. Kaiser 1991; Evrard & Henry 1991). 

The maturing of the ROSAT database has sparked 
much recent interest in testing this result. A large 
cluster sample with a median depth of z ~ 0.1, 
created from the ROSAT All Sky Survey (RASS), 
shows no sign of evolution out to z = 0.3 (Ebeling 
et al. 1997). Castander et al. (1995) presented the 
first look at the high redshift cluster population with 
ROSAT, claiming that the evolution seen by Henry et 
al. (1992) extends to luminosities ~ 10 44 erg s _1 . We 
have recently shown (Collins et al. 1997), using the 
Southern SHARCQ sample of serendipitously detected 
clusters from deep ROSAT PSPC pointings, that the 
number of high redshift clusters is consistent with a 
no evolution model, in direct contrast to Castander 
et al. (1995). Finally, the EMSS sample has been re- 
analysed in the light of new optical and X-ray data, 
which indicates that the evidence for evolution seen 
by Henry et al. (1992) is not statistically significant 
(Nichol et al. 1997). 

In this letter we present the high redshift X-ray lu- 
minosity function (XLF) of the Collins et al. (1997) 
cluster sample. In section || we describe the calcu- 
lation of the XLF from this sample and in section || 
we discuss the results. Throughout this letter we have 
assumed Ho = 50 km/s/Mpc and qo = 0.5 and quoted 
luminosities in the 0.5 to 2.0 keV pass band, unless 
explicitly stated otherwise. 

2. Determination of the XLF 

The cluster sample used in this letter is the South- 
ern SHARC survey, consisting of 16 clusters in the 
redshift range z — 0.3 — 0.7 with a median redshift 
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of z = 0.44, detected in a serendipitous survey of 66 
deep ROSAT PSPC pointings, covering a total search 
area of 17.7 deg 2 . This area is slightly larger than that 
used in Collins et al. (1997) due to the use of a differ- 
ent central source mask. Each pointing satisfies the 
following criteria: exposure time greater than 10 ks, 
|6| > 20°, and S < 20°. Further details of the sample 
are given in Collins et al. (1997) and a description of 
the full survey will be presented in a future paper. 

For each cluster the count rate was measured from 
the background-subtracted count rate image in the 
0.5 to 2.0 keV pass band. The chosen aperture en- 
closes 80% of the light and was determined from a 
convolution of a King surface brightness profile with 
a model of the off-axis point spread function, using 
a cluster core radius of 250 kpc and (3 = 2/3 (Jones 
& Forman 1984). Justification for the adopted values 
of these cluster parameters is given below. After ac- 
counting for the flux falling outside this aperture, we 
used a thermal bremsstrahlung spectrum of temper- 
ature 6 keV (typical of nearby clusters) and a Galac- 
tic absorption model to convert the total count rate 
into a flux, and then into rest frame luminosities, L, 
in both the 0.5 to 2.0 keV and 0.3 to 3.5 keV pass 
bands. We chose to use a thermal bremsstrahlung 
model because cluster X-ray emission is dominated 
by the continuum radiation rather than line emission 
(e.g. Sarazin 1988) and for its computational ease 
of use. To investigate the uncertainty introduced by 
using a single temperature we have used a luminosity- 
temperature relation (Wang & Stocke 1993) to obtain 
iterated estimates of the cluster luminosities. These 
values are within 1% (0.5 - 2.0 keV) and 4% (0.3 - 
3.5 keV) of those obtained above and the difference 
is negligible compared to the Poission errors of the 
photon counts, which has a median value of 10%. 

We have performed extensive simulations to model 
the effect of using an extent criterion to select clus- 
ter candidates. The ability to detect a given cluster 
depends on properties both intrinsic and extrinsic to 
the cluster — the relevant cluster characteristics are 
its luminosity, surface brightness profile and redshift; 
those for the survey are exposure time, background 
count rate and off-axis angle. The simulations as- 
sume a universal cluster X-ray surface brightness pro- 
file, namely the King model with a core radius of 250 
kpc and (3 — 2/3 (Jones & Forman 1984), that is 
independent of redshift. Ideally one would use a dis- 
tribution encompassing the true range of cluster pro- 
files, however the sensitivity and spatial resolution of 
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the current generation of X-ray instruments limit the 
knowledge of such a distribution to low redshift. Our 
choice of profile is consistent with the average prop- 
erties of low redshift cluster samples (e.g. Henry et 
al. 1992; Jones & Forman 1984). As the simulations 
are computationally intensive we limited our analysis 
to a binned representation of the survey with regard 
to exposure time, background count rate and off-axis 
angle. The bin sizes were chosen so as to be small 
enough that the selection function did not vary sig- 
nificantly across a bin whilst being large enough that 
the simulations could be performed in a reasonable 
time. For each bin a range of cluster luminosities 
(0.1 - 3.0 x 10 44 erg s^ 1 ) and redshifts (0.2 - 0.9) 
were used; simulated clusters were generated by ran- 
domly distributing source photons using the King pro- 
file convolved with a model PSPC PSF and back- 
ground photons using a flat background model. The 
number of source and background photons were Pois- 
son distributed about the expected number. 

Multiplying the selection function of a given bin by 
the survey area for the bin and summing over all the 
bins produces il(L, z), the available survey area as a 
function of cluster luminosity and redshift, shown in 
Figure 0. The survey volume for a cluster with lumi- 
nosity L, in a redshift shell z — z\ to z%, is V ulax {L) 1 
defined as 



Knax(^) 



n(L,z) dV{z) dz, 



(1) 



where dV{z) is the volume per square degree at a 
redshift z. Figure || shows equation (|l]) evaluated for 
the redshift shell z — 0.3 — 0.7 using the curves from 
Figure |l| where the solid line indicates the fit used to 
interpolate these points. 

We have used both parametric and non-parametric 
forms to define the cluster XLF; a brief description 
follows and a full account of the fitting procedure can 
be found in Nichol et al. (1997). The non-parametric 
XLF, n(L), was calculated using 



n(L) = J2 



V max (Li) AL 



(2) 



where the sum is performed over all clusters whose 
luminosities, Li, are within ±AL/2 of L. For a model 
XLF <p{L), the parametric form was determined by 
maximising the likelihood L, given by 



J °° V max (L) <f>(L) dL ' 




(3) 



0.4 0.6 
Redshift (z) 

Fig. 1. — The available survey area, as a function of 
redshift, for the Southern SHARC survey. Each curve 
corresponds to a different cluster luminosity, labelled 
in units of 10 44 erg s _1 , and was calculated using the 
simulations described in the text. 

where the product is over the cluster detections and 
Ei(L) is a normal distribution, with mean Li and 
standard deviation equal to the error on Li, which 
models the effect of luminosity errors on the fitting 
procedure. 

In Figures [| and ^ we show the non-parametric 
XLF, for the redshift shell z — 0.3 — 0.7, calculated 
using equations ([!]) and (|J) for both the 0.5 - 2.0 keV 
and 0.3 - 3.5 keV pass bands. Table |l| lists the data 
points shown in these figures. 

We fitted two parametric forms to the data: a 
power law (e.g. Henry et al. 1992), 

d>(L) = KL- a , (4) 

and a Schechter function (e.g. Ebeling et al. 1997), 

(j){L) = K s exp(-L/L m ) L~°" . (5) 

As we are not sensitive to the characteristic luminos- 
ity, L„, of the Schechter function, we used the best 
fit values from Ebeling et al. (1997) — 5.70 x 10 44 
erg s" 1 (0.5-2.0 keV) and 10.7 x 10 44 erg s" 1 (0.3- 
3.5 keV). The likelihood for a given 4>{L) is inde- 
pendent of the normalisation (K or K s ), so we fixed 
this by setting the expected number of clusters equal 
to the number in our sample. This procedure gives 
some weight to luminosities where no clusters were 
detected, producing a steeper slope than the non- 
parametric XLF would suggest (Henry et al. 1992). 
Errors on a and a s were calculated by integrating the 
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Fig. 2. — The symbols indicate the values of V max (L), 
for the redshift shell z = 0.3 — 0.7, calculated using 
equation ([l]) and the curves shown in Figure |l|. Lu- 
minosities are in units of 10 44 erg s _1 and the solid 
curve is the fit to these points. 

normalised likelihood distributions, which are well ap- 
proximated by a Gaussian, out to an enclosed area of 
68.3% and those for the normalisations were found by 
allowing a to take its one-sigma values. The results 
are listed in Table || and shown graphically in Figures 
| and |. 

3. Discussion 

From Figure || it is qualitatively obvious that, for 
L ~ 10 44 erg s"\ our z = 0.3 - 0.7 cluster XLF is 
consistent with that of the low redshift cluster XLF. 
Our best fit faint end slope for the Schechter func- 
tion, a s — 1.77 ± 0.30, agrees with that of the low 
redshift cluster sample of Ebeling et al. (1997), who 
find a s = 1.85 ± 0.09. The power-law slope in the 
EMSS pass band (a = 2.22 ± 0.25) is also consis- 
tent with the slope of the low redshift shell of Henry 
et al. (1992) and the re-worked EMSS sample of 
Nichol et al. (1997), for which a = 2.19 ± 0.21 and 
a = 2.60 ± 0.37 respectively The lack of luminos- 
ity evolution agrees with our recent analysis of the 
redshift distribution of the Southern SHARC clus- 
ters (Collins et al. 1997) and is inconsistent with the 
strong evolution claimed by Castander et al. (1995). 
Our results strengthen, and extend to higher redshift, 
recent claims that cluster properties remain the same 
out to z ~ 0.3 (e.g. Ebeling et al. 1997; Mushotzky 
& Scharf 1997). 

The only other published XLF for X-ray selected 
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Fig. 3.— The XLF of the Southern SHARC in the 
0.5 to 2.0 keV pass band. The solid line is our best 
power-law fit, the dotted line our best fit Schechter 
function and the dashed line is the Schechter function 
fit of Ebeling et al. (1997) to their low redshift RASS 
cluster sample. 

clusters covering our redshift range is that of Henry 
et al. (1992) from 21 EMSS clusters between z = 
0.3 — 0.6, with a median redshift of 0.33. They find 
a much steeper power-law slope of a = 3.27 ± 0.29, 
which they claim is the result of negative evolution. 
Our XLF is consistent with the non-parametric XLF 
of Henry et al. (1992) for luminosities ^ 5 x 10 44 
erg s _1 (0.3 - 3.5 keV), as shown in Figure ||, hence 
any evolution of the cluster population is restricted 
to luminosities J> 3 x 10 44 erg s _1 . Using the Ebel- 
ing et al. (1997) XLF, and assuming 100% detection 
efficiency, we find that the surface density of clusters 
with L > 3 x 10 44 erg s , in the redshift range 0.3 
to 0.7, is 0.083 deg -2 and the expected number of 
such clusters in the Southern SHARC survey is only 
~ 1.5. Therefore the lack of high luminosity clus- 
ters in our sample is not unexpected. In Figure ^, we 
plot equation (6) of Henry et al. (1992) using their 
best-fit values of n = -2.10 and k = 0.029ft. Mpc -1 , 
for z — 0.37, corresponding to the median redshift 
of the joint cluster sample. This analytical model is 
based on the Press-Schechter formalism and is shown 
here simply to provide a physical basis for extrapolat- 
ing the EMSS data and comparing the two samples. 
In the context of general cluster formation theories, 
our result is broadly consistent with the predictions 
of entropy-based models or a low value of the density 
parameter (e.g. Kaiser 1991; Evrard & Henry 1991; 
Henry et al. 1992; Oukbir & Blanchard 1997; Bower 
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Table 1 

The non-parametric XLF for the Southern SHARC survey. 



pass 


band 


log L a 


logn(L) b 




(keV) 










0.5 


- 2.0 


43.50 


-5.52 


(+0.37,- 


-0.45) 


0.5 


- 2.0 


43.80 


-6.26 


(+0.25,- 


-0.28) 


0.5 


- 2.0 


44.10 


-6.61 


(+0.20,- 


-0.22) 


0.5 


-2.0 


44.40 


-7.17 


(+0.25,- 


-0.28) 


0.3 


-3.5 


43.75 


-5.78 


(+0.37,- 


-0.45) 


0.3 


-3.5 


44.05 


-6.51 


(+0.25,- 


-0.28) 


0.3 


-3.5 


44.35 


-6.86 


(+0.20,- 


-0.22) 


0.3 


-3.5 


44.65 


-7.42 


(+0.25,- 


-0.28) 



a L has units of erg s 1 

h n(L) has units of Mpc" 3 AL^ 1 

Note. — Luminosity bins have a constant 
width and were chosen so that each bin con- 
tained at least two clusters. Poisson errors 
(one-sigma) are from Gehrels (1986). 



Table 2 

The parametric XLF for the Southern SHARC survey. 



pass band 


a 


K a 


a s 


if/ 


(keV) 










0.5 - 2.0 


2.15 + 0.23 


2.72 + 0.12 


1.77 + 0.30 


3.64 + 0.16 


0.3 - 3.5 


2.22 + 0.25 


4.99 + 0.77 


1.78 + 0.30 


5.68 + 1.04 



a Units of 10~ 7 Mpc- 3 (10 44 erg s" 1 ) 
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